Searching for vortex structures in high Reynolds number turbulence 
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In experimental study of very high Reynolds number turbulence, we found evidences that 
there are distinguished vortex structures in the intermediate range, that is, between the 
Kolmogorov and Taylor microscales, where they are indeed expected to be present. These 
structures are responsible for the intermittency, and, in the same time, they contribute into 
asymmetry of turbulent statistics, the latter following from the Kolmogorov law. 

PACS number(s): 47.27. Ak, 47.27. Jv 



I. INTRODUCTION 

The presence of organized structures in fully developed 
turbulence was demonstrated recently both in physical 
experiments, and in numerical simulations. Most of the 
studies are devoted to vortex tubes (or "worms", or" 
sinews"), see, e.g., (|]- @- 

Generally, these organized structures may be consid- 
ered at the same time as a source (or manifestation) 
of intermittency, because normally they correspond to 
a small dimension (of the size of Taylor microscale, or 
Kolmogorov microscale), although they are stretched in 
the other directions over a "macroscopic" length. 

One remarkable feature of these vortex structures is 
that they present not only the intermittency of turbu- 
lence, but also contribute into the asymmetry of the 
PDF's M, |Q. As was suggested in 0, the asymmet- 
ric statistics of the turbulence is indeed related to the 
intermittency. This connection was summarized in the 
so-called ramp-model, the latter being only empirical. 
The considerations of the formation of vortex structures, 
when the intermittency and asymmetry appear simulta- 
neously, have resulted in the development of the ramp- 
model into a more sophisticated version which is now 
called the bump-model. 

The data used in this paper are based on Taylor's hy- 
pothesis, and therefore do not provide direct measure- 
ments of vortices. However, these vortices should be 
manifested in different ways, such as appearance of in- 
termittency, etc. We provide here evidences that these 
vortices are present in fully developed turbulence. We be- 
lieve that the disadvantages of necessary indirect studies 
of vortices presented in this paper, are outweighed by 
the fact that we are studying a very high Reynolds tur- 
bulence. 

The paper is organized as follows. Section [n]is still in- 
troductory, presenting (mostly known) notions of turbu- 
lence statistics, related to the vortices. The main purpose 
of this section is to explain, what can we expect to ob- 
serve from experimental data if the vortices are present, 
and also to introduce denotations we will use in the rest 
of the paper; and thus some repetition of previous results 
seems to be inevitable. On the other hand, the bump- 



model is introduced in Sec. [II 13. Section [II starts with 



suggestions how to improve the two extreme cases, or 
models, considered in Sec. [fl[ The answer is in the mid- 
dle: an intermediate model is suggested. The rest of the 
section is devoted to comparison of this model with the 
experimental data, namely, with the observed structure 
functions of low orders. Section ^ presents evidences 
of existing vortices. Namely, it is shown that the inter- 
mittency is substantially enhanced in the intermediate 
range: from studying the fl atness of the velocity incre- 
ment s stat istics, Sec. IV A , and from the box counting, 
Sec. [IV B| . An extensive discussion of both extreme case 
models, and also of the intermediate model is given in 
Sec. 0. It is argued that this model may give an insight 
into the question of the intermittency formation, which 
proves to be intimately connected with the asymmetry 
of turbulence. Finally, concluding remarks are presented 
in Conclusion, Sec. Nl. 



II. BASIC IDEAS AND DENOTATION 



We consider structure functions, 
S n (r) = {(A r u) n ), 



(1) 



where A r u = u(x + r) — u(x), and n is an integer, and 
generalized structure functions, 



S q (r) = (|A rU |«), 



(2) 



where q > 0. 

We mention an exact result, obtained directly from 
Navier-Stokes equation, 



S 3 ( r ) = -Lr + 6u^-S 2 (r), 
5 or 



(3) 



[gj, and in inertial range where the viscosity term can 
be neglected, we recover the 4/5-Kolmogorov law, 



5 3 (r) 



4 

-—sr. 
5 



(4) 
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Data are from atmospheric turbulence measurements, 
about 35 meters above the ground, obtained by hot- 
wire anemometer, mean wind speed was 7.6 m/sec, root- 
mean-square velocity 1.3 m/sec. Using Taylor's hypoth- 
esis and local isotropy assumption to obtain dissipation, 
one obtains the Taylor microscale Reynolds number to be 
9540 and the Kolmogorov microscale r\ to be 0.57 mm. 
The data were sampled at 5 kHz and the file consisted 
of 10 million data points (courtesy of Sreenivasan) . The 
data are treated in the spirit of the Taylor hypothesis, 
that is, the time series is treated as one-dimensional cut 
of the process (for more detail, see 



A. No intermittency. Asymmetry in K41 

Consider classical picture of turbulence due to K41 
jL2[. That is, the turbulence is presented with an en- 
semble of multiscale cells, from integral scale £ to small- 
est Kolmogorov scale rj ~ £/i? 3 / 4 , where R w vl/v, the 



Reynolds number, and v = y (v' x + Vy + vf ) f the mean 

square velocity, with |A r tt| ~ r 1 / 3 . Then, in inertial 
range, i.e., at r\ < r < £, 



S q (r) 



r 9/3 



We define the Taylor microscale 
A 

see, e.g., Then, 

I 



A 



R 1 / 2 

It follows from (pf) that 




and R\ = 



A = 1&!%RV 2 



x i 



(5) 



(6) 



(7) 



(8) 



see, e.g., Q, formula (3.2.18), we use this expression 
below. 

Let us proceed to the asymmetry, dictated by the Kol- 
mogorov law (|J). Denote Ai ; 2,3 - three eigen values of 
the symmetric part of the matrix djVi, and suppose that 



Ai < A2 < A3. 



(9) 



It is known Jl5j, that the quantity defining skewness of 
d x u field, namely ((d x u) 3 ), is = (8/35)(AiA2A3), and 
therefore we calculate (A1A2A3). As Ai + A 2 + A3 = 
(due to the incompressibility), (0) can be replaced by 



y < A 2 < A 3 . 



(10) 



We thus consider conditional mean value (A3 is fixed) 
(AiA 2 A 3 |A 3 ) = -<(A 2 + A 3 )A 2 A 3 ) = 



-A3/2 



p(A 2 |A 3 )[-(A 2 + A 3 )A 2 A 3 ]<iA2 



(11) 



where p(A2|A 3 ) is conditional probability distribution of 
A2. The quantity A1A2A3 as a function of A2 is highly 
asymmetric and mostly negative in the range (0). If, 
for example, we suppose that p(\2\^3) is an uniform dis- 
tribution, i.e., 

1 

p(A 2 |A 3 ) = 



(3/2)V 

which at first sight seems to be quite a realistic assump- 
tion, then, according to 0), (A1A2A3IA3) = -Af/2 < 0. 
In fact, the contribution of the negative part of A1A2A3 
is 10 times larger than the positive part; that is, in or- 
der to make the integral in (11) positive, the distribution 
^(^2^3) should be highly asymmetric. Namely, p(\2\^3) 
should be more than 10 times larger in the negative range 
of (|T^) than in the positive. As this kind of assumption 
concerning the PDF p(X2\Xs) does not seem to be reason- 
able, it follows from these considerations that ((d x u) 3 ), 
and therefore the skewness should be negative. 

In a more natural way, one may consider the strain 
parameter, 



A = 



Ai — A2 _ 2A2 + A3 
Ai + A 2 ~ 



A a 



introduced by Moffatt et al. and, according to (|lC|), 
< A < 3. Then, 

A1A2A3 = ^(1-A 2 ). 

The interval where this expression is positive (i.e., < 
A < 1) is twice as short as where it is negative (i.e., 1 < 
A < 3), and therefore the mean value may be expected to 
be negative. And indeed, for an uniform PDF p(X\Xs) = 
1/3, we recover 

(AiA 2 A 3 |A 3 ) = ^ [ p(X\X 3 ){l - X 2 )dX = -X 3 /2 < 0. 
4 Jo 

These considerations appear to be "too general", and 
therefore suspicious, because they do not incorporate the 
equation of motion. And indeed, these considerations are 
bias. Recall that A3 (which is > 0) was previously fixed; 
if we now fix Ai(< 0), then 



Al < A 2 < -y, 



(12) 



cf, (|To|). Calculating now (A1A2A3IA1) = — (AiA2(Ai + 
A2)), analogously to (|ll|), we obtain (A1A2A3IA1) = 
-A 3 /2 > 0! 

Finally, if we no longer require (^|), i.e., if we do not re- 
quire any ordering for the eigen values, then (A1A2A3) = 
0. Indeed, all probability distributions are expected to 



2 



be symmetric a priori in respect to the transformation 
^1,2,3 — > — Ai,2,3, and therefore (A1A2A3} = 0. Thus, 
these "general considerations" fail to account for the 
asymmetry, because they do not incorporate the equation 
of motion. In the framework of K41, the asymmetry is 
needed to account for the energy transfer to small eddies 
]13j , the asymmetry by itself being expressed through the 
Kolmogorov law (Q). However, the theory does not in- 
dicate what dynamical processes lead to the asymmetry. 
Another setback in this picture is that the turbulence 
is completely non-intermittent, in contradiction with ex- 
perimental data. 

In particular, the asymmetry dictated by the Kol- 
mogorov law, and observed (negative) skewness, are not 
supposed to be related to the intermittency. That is, the 
PDF for A r u is asymmetric, however, the asymmetry 
has nothing to do with the tails of the PDF, that is with 
possible intermittency. In other words, the asymmetry of 
the PDF is expected to be manifested in the core. This 
"ideal" PDF is constructed in O], in order to compare 
with observed properties of turbulence. The PDF, 



p(u, r) — T{u), u 
has the following properties, 



A, 



((A.u) 2 ) 1 ^ 



(13) 



I(u)du = 1, / ul(u)du = 0, / u 2 X(u)du = 1 



and 



u 3 X(u)du = — 



4 1 



5 C 3/2 ' 



(14) 



where C2 is the Kolmogorov constant. But, most impor- 
tantly, the PDF is constructed from two Gaussian func- 
tions, so that it does not contain any tails. We will use 
this "ideal" PDF below; we note for now that the exper- 
imental asymmetry is related to the intermittency. In 
particular, the tails of real PDF giv e a substantial con- 
tribution to the Kolmogorov law |L1|] . 

In addition, the experimental PDF p(u, r) is not self- 
similar, that is, it is also a function of r (as well as of 
u, cf. (|l3|)), see, e.g., Q. This in effect implies presence 
of intermittency. Direct study of intermittency can be 
provided by cumulative moments, 



So,c( r ) = / p{u,r)du, 



(15) 



and 



S-, c {r) = J p(u,r)du, S +:C {r) = J p(u,r)dii, (16) 
and by the tail moments, 



So,t(r) = I p(u,r)du + / p(u, r)du = 1 - S , c=t (r), 
) Jt 



and 

/ — t roc 
p(u,r)du, S + . t (r) = p(u,r). (18) 
-00 J t 



(we will consider c = t = 3, 



If c and t are not small 
and c = t = 4), then (15) presents major events pro- 
duced by the core of the PDF, and So,c should be close 
to unity, or rather Gaussian value with which it will be 
compared. We denote the cumulative and tail moments 
for Gaussian and "ideal" distributions with letters G and 
X correspondingly. For example, So tC (r) for Gaussian dis- 
tribution is Go jC) and SL^r) for ideal distribution is 2L ,{, 
etc. And, of course, So, c should be independent of r, in 
case of self-similarity. On the contrary, substantial devi- 
ation from unity, and, what is more important, depen- 
dence on the distance - increasing with growing r, would 
indicate intermittency. In fact, super-Gaussian deviation 
from unity suggests that the tails give considerable con- 
tribution to the PDF, and this contribution is decre asing 
with growing distance. We will see in Sec. IV B that 
that is what is indeed observed. 

As to S± jC (r), they represent a more subtle way to 
measure intermittency, and, of course, asymmetry. These 
quantities are expected to be close to (and smaller than) 
1/2, and 5+ jC being slightly larger than 5_ iC . 



B. Extreme intermittency 

There is a geometrical interpretation of (A1A2A3} < 0. 
It implies that typically Ai < 0, while A2,3 > 0, that 
is, the inflow is in one direction, 1, while the outflow 
proceeds in two directions, 2,3. Statistically, these flows 
would generate vortex sheets rather than vortex ropes 
|l6| . However, it became clear that the presence of the 
vortex itself changes the principal axes and eigen values 
|&B. Consider an axisymmetric cell of a scale £. At the 
axis, there is an inflow in two directions, 1,2, that is, 
to the axis, and outflow in one direction, 3, parallel to 
the axis. Therefore, locally A1A2A3 > (Ai,2 < and 
A3 > 0). To be more specific, in cylindrical coordinates 



v = { - 



0, 2az}, 



(19) 



a = — Ai = — A2 = A3/2. Suppose now that there is a 
symmetric vortex present, with initial scale £. Then the 
vortex is stretched by the strain motion (|l^), according 
to equation 



dtUJ z — ar'd r >uj z = 2au> z -\ — -d r 'r' d r 'Uj z . (20) 



This stretching proceeds until viscous effects come into 
play, resulting in Burgers vortex, 



(17) 



u z (r') = w z (0)e 



-r' 2 /A 2 



(21) 
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FIG. 1. Schematic presentation of bump model, and its comparison with the ramp-model, (a) Presentation of d x u for the 
strain (^), solid line. The vortex is superimposed on this plot (dashed line), (b) Resulting motion, the strain plus the vortex, 
(c) More realistic presentation of the same as in (b). (d) Corresponding motion u. Note typical bump. The dashed-dotted 
line presents a bump in a "wrong place", that is, this situation is statistically unlikely, (e) d x u for the ramp- model, and (f), 
corresponding motion, which is indeed ramp-like. The dashed-dotted line depicts another ramp, which is supposed to appear 
statistically rarely. 



4 



and 



where 



V<f> 



_ ^(0)A 2 / 
2r' \ 



V a 



(22) 



(23) 



As a ~ i^/A A ~ £i? -1 / 2 , and therefore A coincides 
with Taylor microscale (^) . In order to estimate lu z (0) we 
note that, neglecting viscosity, equation (|2^) conserves 
the quantity 



,r'dr' 



and hence we estimate, lo z ~ 1/f*' 2 , that is, 



Ws(0) « = ^-R) 



(24) 



(25) 



where u>i is initial vorticity. 

It was shown in that the combination of the strain 
motion ( [l9| ) and the vortex results in negative skewness. 
In order to get some simple estimates, we present d x u 
corresponding to this motion in Fig. 0. Figure 0(a) 
presents the strain (|l9| ) with a = 1. The probability 
p(d x u > 0) = 1/3 (where d x u = 2), while p(d x u < 0) = 
2/3 (and d x u — —1), and therefore 



(d x u st ) = 2p(d x u > 0) - lp(d x u < 0) = 0, 



and 



Figure 0(c) depicts essentially the same as Fig. |1|( 
but it suggests a more realistic presentation of d x u. The 
integral of this latter function, that is, the velocity u itself 
is depicted in Fig. 0(d). Note characteristic bump which 
appears due to the vortex. Panels (e) and (f) are to be 
compared with (c) and (d); they present d x u and u for 
the ram p mo del ( the d etails of this comparison are given 
is Sees. EVBl and[VA|). 

Imagine now a statistical ensemble of the cells of the 
scale £, so that the ensemble as a whole is statistically 
isotropic. Suppose that each of the cells contains a Burg- 
ers vortex ( pl| ) with u>i « a. Then, 

({d x uf)~^R^~u}R~(^) 2 R : (28) 



so that 



v((d x uf)~^~e, 



(29) 



as in K41 theory. 

The Kolmogorov law (0) is also valid, of course. How- 
ever, the "inertial range" for this ensemble, starting at 
the scale £ is cut off already at the Taylor microscale A, 
rather then Kolmogorov microscale rj, because the linear 
in r term in (Q) is larger than the viscous term only at 
r > A. At r < A 5a (r) ~ r 3 , reaching at r = A the am- 
plitude 53(A) ~ eX, as should be according to (0). In ad- 
dition, the skewness, ~ 53(A)/ A 3 ~ — e/A 2 w -{vi/£) 3 R 7 
which coincides with the estimate fl27|). 

So far, this ensemble looks similar to the K41 the- 
ory. However the intcrmittency in this model is much 
too high. For example, {(d x u) 4 ) ~ lu 4 R 4 X 2 /£ 2 ~ u 4 R 3 , 
so that the flatness, 



((d x u st ) 3 ) = 8 P (d x u > 0)-lp(d x u < 0) = 8/3-2/3 = 2 > 0. 

That is, the skewness is positive for the strain motion, 
as mentioned above. The vortex corresponds to a sharp 
odd function concentrated on a small scale A. Therefore, 
{d x u\) = and {{d x u\) 3 ) = 0. Nevertheless, the sum 
of these motions, depicted in Fig. 0(b), may possess 
negative skewness. Indeed, 

(d x u st + d x u x ) = (d x u st ) + (d x u x ) = 0, 



F(0) 



((d x u) 4 ) 
{(d x uW 



R, 



(30) 



which is many orders of magnit ude ab ove the experimen- 
tal value (see, e.g., below Sec. [V A). 



III. 



COMPARISON WITH EXPERIMENTAL 
DATA 



while 



((d x u st + d x u x ) 3 ) = ((d x u st ) 3 ) - 3 • l((d x u x ) 2 ) (26) 

(cf. Q), which becomes negative if the second term on 
the right-hand-side prevails. In particular, for the strain 
@ plus vortex @, 



{(d x u st + d x u x ) 3 ) = 2a 3 ~ 3aLo 2 e R 2 -^ « -3a^ 2 i?, (27) 

where we used the estimate (p^), and the fact that the 
vortex occupies only X 2 /£ 2 part of the space (naturally, 
we also suppose that tug is not small compared with a). 



Further discussion of these two models and the ramp 
model is provided below in Sec. |y|. It is clear for now 



that both extreme cases: no intermittency 
treme intermittency 



[IB 



II A 



and ex- 



are not realistic, and it should 
be accepted something in between. Let us accept the 
K41 theory as a "first approximation" , and we will try 
to modify it by introducing finite intcr mitte ncy. Then, in 
addition to what is described in Sec. II A. the cells are 



generating vortices as in Sec. [IB above. The question 



is, what are the initial magnitudes of Lo r l It is clear that, 
for example, u>i should be small compared with vg/£: oth- 
erwise the intermittency is too strong. Another possibil- 
ity is to suggest that uj£ ~ vi/£, however, only a small 



5 



fraction of the cells are generating these vortices. In addi- 
tion, there is nothing special in the large scale eddies; in 
fact, the cells with scales r < I are generating the vortex 
ropes as well. Namely, the scales of generated vortices 
by the cells of the scale r are, 



I 2v 

v r /r 



= A 



1/3 



(31) 



cf. (|23J). Each cell is able to generate a vortex rope as 
long as S r < r, i.e., according to (^T|), when r > £/R 3 ^ ~ 
•q, that is, in the whole inertial range. Thus, the vortices 
are generated in each cell, however, the scales of these 
vortices, S r , arc, 



i] < S r < X. 



(32) 



We can call this range intermediate , and it is in this 
range where the vortices are indeed observed H . 



We see that now the vortex is aligned along the inter- 
mediate eigenvector. This feature was initially suggested 
and shown in ||, and confirmed both in numerical sim- 
ulations ID , and in laboratory experiments | fl7| . 

In addition, the eigen values ( J33] ) result in correct sign 
of the skewness (i.e., correct asymmetry of the PDF). 
Indeed, 



A1A2A3 = 2a 



dr' 



r' 



<0, 



if the amplitude of the vortex is large enough. This 
formula is consistent with the estimate (|27|). Indeed, 
dv^/dr' — v^/r' ss u>iR 2 , and averaging this expression 
over the space would result in a factor ~ 1/-R, because 
the vortex is occupying only 1/R fraction of space. As a 
result, we recover (p7|). 



A. The third moment 




)K )K X experimental data 
fit 

-A, w (r) 

-A 3 w (r) 



In order to understand what impact has this model 
on the structure functions, we start with the third or- 
der structure function. As mentioned, for the extreme 

^3 



[IB, S; 



II A 



r" for r < A, while 
, 53 ~ r 3 is valid 
incd above in this 



intermittency case, Sec. 
for no-intermittency case, Sec. 
for r < rj. Naturally, the picture out 
section suggests intermediate situation. 

Figure |2| depicts experimental Ss(r), and it is com- 
pared with wh at we would expect for the ensemble de- 
scribed in Sec. II B, namely, with the following function 



constructed for this purpose 



r/\ 



r/r, 



FIG. 2. Experimental and theoretical (somewhat simplified 
and sketchy) third moment, Ss(r). 



We note that, according to Sec. p B| , the vortex is 
stretched along the stretching axis (with maximum value 
of Ai 2 3)) of course, that is, parallel to the z-axis. How- 



ever, as mentioned at the beginning of Sec. II B, when 



the vortex is formed, the eigen values are changed. While 
A3 remains the same = 2a, the other two eigen values are 
now 



Al,2 



-«4 



~dV 



V<j> 



where is defined in (|22]), see, e.g., 0. We now have to 
arrange the eigen values in order (H) to obtain 



Ai = 



-a 

2 



dv,p 
~dV 



v 



A 2 = 2a, 



-— retanh . 
5 V A 



The experimental £3 (r) is also compared with 

2 



4 , r 
-—re tanh — 

5 \ V 



which is what one would expect from K41. Note that 
all three curves, one experimental and two constructed, 
match at inertial range, i.e., at r > A. In order to do this, 
the constructed functions are matching experimental lin- 
ear fit with exponent £3 = 0.99 ± 0.01, quite close to the 
Kolmogorov law (for which £3 = 1). The plots are given 
in Kolmogorov microscale, that is, in r/r], and there is 
also r/A scale in the figure for comparison. It can be seen 
that the inertial range indeed starts at r ~ A, and thus 
the intermediate range (|3^) is clearly visible. Besides, the 
experimental plot is situated between no intermittency 
moment, 1S3 \r), and extreme intermittency, S^\r), - 
as expected. 



A3 = 



dv<f, _ U0 
dr' r' 



(33) 
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FIG. 3. Panels (a) through (d): the generalized structure functions: starting from the first and ending with the fourth. 
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B. Other moments 

The inertial range is associated with scaling, that is, 
in the log-log plotting the moments are presented with 
straight lines. Figure || illustrates generalized structure 
functions of different orders, the first, the second, the 
third and the fourth (and, of course, the structure func- 
tions of even orders coincide with the generalized struc- 
ture functions). The intermediate range ( |32|) is easily 
separated from the inertial range (where the plots be- 
come straight lines). 

We may interpret this separation of the intermediate 
range by the fact that this range contains vortices. This 
interpretation is not unequivocal however. It is also pos- 
sible that the depletion of the structure functions in this 
range, as compared with the K41, is due to the action 
of viscosity. Indeed, the viscosity is still acting even at 
r 3> 77, in spite of the fact that the local Reynolds number 
is there ^> 1. The fact that the viscous effect is stretched 
up to r ~ A might be purely coincidental. 



stormy) events, the odd order structure functions are sub- 
ject to substantial fluctuations. 

Essentially the same conclusions can be drawn from an- 
alyzing Fig. ^], where the dimensionless ratios of different 
moments to the second order structure function are de- 
picted. It is clear that the intermittency is substantially 
increased in the intermediate range. In particular, for 
the ratio S\{r) / S2{r) 1 ^ 2 , the presence of intermittency 
implies that the ratio is below the Gaussian value, see 
Fig. |^(a), rather than otherwise for higher moments; 
and still, the intermediate range is clearly distinguished 
by enhanced intermittency. It is also clear from panel 
(b) that the skewness substantially exceeds that given 
by the "ideal" PDF only in the intermediate range. In 
the inertial range, the skewness is decreasing systemati- 
cally with increasing distance r, although the deviation 
from "ideal" behavior is not conspicuous (and the scaling, 
if any, is rather poor). Finally, maximal flatness (panel 
(d)) reaches the value 32.28, well below what would be 
expec ted f rom extreme intermittency model described in 
Sec. ETb. 



IV. INTERMITTENCY 



B. Box counting 



A. Flatness 

A real test to check if the intermediate range is related 
to vortices is achieved by the measuring the intermit- 
tency. Traditionally, the later is estimated by measuring 
the flatness, see, e.g., More specifically, the flat- 

ness, and related quantities, like generalized flatness, are 
compared with corresponding Gaussian values. Figure || 
presents dimensionless ratios of different moments to the 
first moment. They can be called generalized flatness. 
We see that, first, the intermediate range is clearly dis- 
tinguished from inertial, the intermittency growing sub- 
stantially in the intermediate range. The deviation from 
Gaussian and "ideal" behavior is especially dramatic in 
Fig. ^(d) , which is hardly surprising taking into account 
that it represents a relatively high moment. 

Second, although the deviation from Gaussian and 
"ideal" values in inertial range is not big, (a factor 
of 2 or 3), there is a noticeable scaling for about two 
decades. Recall that in the framework of K41, accord- 
ing to (||), all these ratios would be flat, that is, e.g., 
/ Si(r) 2 = const, etc. If however the structure 
functions have anomalous scaling, S q (r) ~ r^ 3 , where 
C q ^ q/3, then these ratios would possess some scaling, 
which is indeed observed. 

Finally, what may be called generalized skewness, de- 
picted in Fig. ^(b), is presented by rather scattered 
curve. This might be explained as follow. The third mo- 
ment of the structure function, appearing only due to 
the asymmetry (as opposed to the third order general- 
ized structure function, depicted in panel (c)), contains 
quite a substantial contribution coming from the PDF 
tails, see, e.g., Eoi. As the tails correspond to rare (but 



Zeroth moment of structure function, that is box 
counting, are very useful to study because, unlike high 
moments, the statistic is as good as one can get, although 
quite often the zeroth moment provides only trivial re- 
sults. In order to extract some useful data in studying 
intermittency, one can either analyze the cumulative (and 
tail) moments, or separate negative values from the pos- 
itive, and to study them separately (as, e.g., in |ll|). 

Panel (a) and (b) in Fig. [6] correspond to experi- 
mental cumulative and tail moments, defined in (|l^) and 
( |l6| ) correspondingly, and we compare them with the- 
oretical values. As mentioned in Sec. II A, for self- 
similar PDF's, the cumulative moments are independent 
of distance r, and they are expected not to deviate much 
neither from Gaussian nor from "ideal" asymmetric dis- 
tribution, - these two happen to be close to each other, 
as seen from panel (a). On the other hand, the contri- 
bution of the tails, that is, intermittency, if present, is 
increasing with decreasing distance. That implies that 
the cumulative moments, that is, the cores of the PDF's 
give a smaller contribution for decreasing distances. Fig- 
ure ^(a) shows that this is indeed the case. It can be 
seen that, first, the cumulative moments noticeably devi- 
ate from both Gaussian and "ideal" values. Second, this 
deviation is higher in the intermediate range, and that 
is the way this range is easily separated from the iner- 
tial range. One can see analogous trends from panel (b), 
where tail moments are depicted. The tails definitely 
give larger contribution for small distances, and this is 
especially true for the intermediate range, which is again 
noticeable. Note also that the tail parts for t = 4 are 
substantially larger than Gaussian values, especially for 
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FIG. 6. Cumulative and tail moments of zeroth order, (a) The probability of |«| < 3, and of \u\ < 4 for different r's, which 
coincide with So,c=3(r) and <So,c=4(r) correspondingly. These quantities are compared with correspondent Gaussian and "ideal" 
values, Go, c =4 being almost indistinguishable from unity, (b) Tail moments, with respect to correspondent Gaussian values, 
and for the same values t = 3 and t = 4, are depicted, (c) The probability of u < 0, corresponding to S- ;C ^-oo(r), probability 
of — 3 < u < 0, or S-,c=3{r), and of — 4 < u < 0, or S-, c =4,{r) are depicted, and compared with corresponding Gaussian values. 
This time, G_ jC =4 is practically indistinguishable from 1/2. (d) The probability of u < — 3 (or 5_,t = 3(r)), and of u < — 4 (or 
S_ j t=4(r')) are depicted in respect to their Gaussian values. 
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the intermediate range. 

We now proceed to the analysis of the ± moments. 
It is clear that there are several common features for 
ramp model and bump model, that is, common features 
in p anels (c) and (e) of Fig. [I]. As mentioned in Sec. 



II B| , for both cases, (d x u) = and ((d x u) 3 ) < 0. This 



actually implies that the positive part of the distribution 
occupies larger space than the negative part. This is seen 
directly from panel (e), and can be calculated for the 
function depicted in panel (c), i.e. calculated directly 
from the function d x u st + d x u\ defined in Sec. [II B| . In 
other words, S-^^ooir) should be less than S+, c ^oo(r). 
As, on the other hand, S-^^ooM + S+.c^ooM = 1, by 
definition, we have S'-.c^ooM < 1/2, and this quantity is 
typically increasing with growing r, while S+, c ~*oo(r) > 
1/2, decreasing with r. This is exactly what is observed 
in panel (c) of Fig. ^, although the curves are somewhat 
scattered. In spite of the latter, there is no doubt that 
S- tC (r) < S + ^ c (r), for all distances and all values of c. 

The box counting <S±,c=3(f) and 5± jC= 4(r) behaves 
analogously. The intermediate range can be again distin- 
guished in all of the curves depicted in panel (c), although 
the range is not that pronounced. This range can be 
better noticed from the tail moments, depicted in panel 
(d). It is clearly seen that the tails are stronger at small 
distances, and the intermediate range clearly noticeable. 
Besides, the negative distribution is definitely more sin- 
gular than the positive, or SV, 4=3,4(7") < S— ,t=3,4(r) : and 
that is what is seen from this panel. This difference in 
singularity strength is obvious from the ramp model, Fig. 
|l|(e), where the positive distribution is not singular at 
all. It is also seen from Fig. 0(c), that is, for the bump 
model, that the negative peak reaches higher absolute 
values than the positive peak, and therefore the nega- 
tive singularity is expected to possess higher strength. 
Finally, again as in panel (b), the tails deviate quite sub- 
stantially from Gaussian values for t = 4. 



V. DISCUSSION OF THE MODEL 



A. Intermittency formation 



Let us return to the stretching process described in Sec. 



II A . In more general case of axisymmetric stretching we 



have, 
~dt 



+ 



1 dr' 
r' dr' 



d_ 

' dz 



V4, 



dz 2 



d_}_d_ 

dr' r' dr' 



V0. 



(34) 



In the vicinity of stagnation line r' = 0, this equation is 
reduced to ( |20| ) if operator Vx acts on (|4|). This equa- 
tion is linear in respect to u^, and, if v r < and v z are given, 
then the vortex stretching, described by equation (|34|) 
can be considered kinematic. In fact, this equation coin- 
cides with that for A^, </>-component of magnetic vector- 
potential, i.e., B r i = —dzAcf), B z = r'~ 1 d r >r' A$, and B is 



magnetic field strength. The induction equation for A^, 
coinciding with (|34|), is kinematic if the velocity is given. 

The equations for r and z-components of the velocity 
field read, 



dv r 
~di 

1 dp 
p dr' 



+ 



d_ 

dr' 



d 



r' 



d 2 d_]_d_ ,' 

dz 2 dr' r' dr' 



v r ,+F r ,, (35) 



and 



dv z 
dt 



+ 



d_ d_ 

r dr' z dz 



1 dp 
p dz 



+ 



dr' r' dr' 



(36) 



where p is density, and F is external forcing. 

In the vicinity of stagnation line r' — > where there 
is the strain (|19|) motion plus the vortex, all the terms 
in equation (|35|) are functions of r' only, and therefore 
they can be compensated by the pressure. Analogously, 
all the terms in ( |36] ) are functions of z only, and again 
are compensated by the pressure. This means the "the 
cell", i.e., the r and z-components of the velocity can be 
considered as given (and steady). In general, however, 
V& depends also on z, and therefore the third term in 
( p5| ) gives nontrivial contribution to the equation. Nev- 
ertheless, as mentioned at the beginning of Sec. Ill, the 



velocity corresponding to the vortex is supposed to be 
small compared with the strain motion. In other words, 
the contribution of this term can be neglected. Besides, 
the external forcing is always able to make the cell steady. 
In any rate, if this term is neglected, then these two equa- 
tions (p5]-p6[) are independent from existing vortex, and, 
from the point of view of the vortex generation, the iy 
and ^-components can be considered as given. 

In this - somewhat simplified - picture, the vortex gen- 
eration responsible for appearance of the intermittency, 
is kinematic, that is, passive. We may expect therefore 
that there is some analogy between this generation of in- 
termittency and that of a passive scalar gradient. It is 
known that the latter is shown to be intermittent even if 



the background, that is, the given velocity is not |18|. 

Generally, we have the following energy balance equa- 
tion, 



dt 



= -W(V x v) 2 + (Vx v w ) 2 ), 



(37) 



where u w stands for the generated vortex. The first 
term on the right-hand-side corresponds to regular Kol- 
mogorov cascade, and corresponding energy dissipation, 
and the second one - to the dissipa tion o f the vortices. In 
the no intermittency model, Sec. 



II A, the second term 



is absent. O n th e contrary, in the extreme intermittency 
model Sec. II B, the first term can be neglected, because 
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of the large size of the cells, while the second makes up 
for the Kolmogorov dissipation if the initial amplitude of 
the vortex ui£ is comparable with ott, see (p9|). In a more 
realistic model described at the beginning of Sec. H , 
both dissipation terms contribute. More specifically, the 
Kolmogorov cascade prevails, and the generated vortices 
give a small contribution to the dissipation (otherwise 
the intermittency is too strong). Let lu£ ~ fote, where 
/ < 1, then 



(38) 



cf. (28). It is because of the smallness of / 2 that the 
process of vortex generation can be considered kinematic. 

Presumably, the most vulnerable part of this picture is 
that the parameter / is essentially arbitrary. We can sug- 
gest some speculation about the origin of this parameter. 
Namely, the process of vortex formation, and vortex ex- 
istence, are transitional. Indeed, the vortex cannot exist 
longer than the cell itself, where the vortex is generated. 
According to K41, a cell of the size I exists only for a 
time ~ i/vi. On the other hand, 1/ai ~ £/vg as well. 
In other words, the vortex generation time, 1/ai, and 
the life time of the cell, and therefore the life-time of the 
vortex itself, are comparable. Of course, these are just 
rough estimations. There is some dispersion both in the 
life-times, and in the generation times. We may expect, 
that a cell should be persistent, in order to have time 
enough to generate a vortex. That is to say that only a 
fraction of cells (which are persistent) would give rise to 
vortices. Another restriction arises from geometry. We 
may expect that there are vortex per turba tions in each 
cell. However, as described in Sec. IIB, a cell would 



generate a vortex only if the initial perturbation is ax- 
isymmetric, with the axis parallel to the stagnation line 
of the cell. Generally, a perturbation, even with u>e ~ a , 
would not be exactly the way described in Sec. IIB. 



Namely, it would not be expected to be axisymmetric at 
all. Considering again the interaction to be kinematic, 
we may expect that only axisymmetric part of the initial 
perturbation, - being generally a small fraction of the 
perturbation as a whole , - is transformed into a vortex, 
considered in Sec. II A. Another solution to this prob- 



lem can be seen in the fact that the generation of the 
vortices does not necessary proceeds in the largest cells. 
Of course, flatness according to (^0|) ~ R is much too 
high. If, however, the generation proceeds from a cell of 
a size r, then, the flatness would be 



F(0) ~R r = R (J) 



4/3 



< R. 



Another difficulty in this model is that it predicts the 
presence of intermittency only in the intermediate range, 
(E32|) in Sec. 



III. However, we can see from practically all 



the plots that the intermittency is present in the inertial 
range as well, although, being less intense. It can be 
seen, for example, from Figs. || and |[ where flatness 



is depicted, that the latter has a decent scaling in the 
inertial range (for two decades!). Nontrivial behavior in 
inertial range can be also observed in Fig. EL especially 
in panels (a) and (b). The tail events are particularly 
noticeable for the tails with t = 4, and their deviation 
from Gaussian is substantial in numbers, as seen from 
panel (b) , although these events are less pronounced than 
in intermediate range (and the curve is more scattered). 

We may suggest the following explanations. First, as 
noted above in this section, the vortices are transitional 
events, being generated at large scales, and then, with 
diminishing scales they reach their final scale A. Thus, 
at some point in time, each vortex has a size comparable 
with r, and therefore it may contribute to the distribu- 
tion. Second, strictly speaking, the structures of smaller 
(and larger) than r scales do contribute to the struc- 
ture functions taken at a distance r. Therefore one may 
expect that the structure functions for distances r corre- 
sponding to the inertial range would reflect the intermit- 
tency of smaller scales, or possibly the intermittency of 
the intermediate range only. 



B. Ramp model versus bump model 

As we noticed in previous section, there are several 
common features between these two models. Here, we 
focus on differences. First, the estimates (p?j-p7|) differ 
for the ramp model. The first moment is still zero, of 
course, (d x u ramp ) = u> • (r - S) - [(r - 5)/S\lo -5 = 0. 
Here S is the size occupied by the negative part, and 
therefore r — 5 corresponds to the remaining positive part, 
the amplitude of the positive part being w, see Fig. |l|(e). 
The third moment is negative. Indeed, 



uj) S < (39) 



{(d x U r a7np) 3 ) = W 3 (r - S) 



which is negative if d < r. Therefore, the skewness is 
negative as in (|2^ - p7|) , but the parameters entering this 
expression are different. 

The main conclusion about the ramp model is that the 
strength of the negative singularity is higher than that of 
the positive, 



(40) 



where are generalized dimensions, JtJ]. In the simpli- 
fied picture illustrated in Fig. |](e), the positive part is 
not singular at all, that is, D+ = 1. On the other hand, 
the bump model does suggest that the positive part is 
singular, although the inequality ( pio| ) is satisfied: it is 
seen from Fig. [j](c), and follows from discussion at the 
end of previous section. The singularity of positive tails 
is indeed observed, see Fig. ||(d), and that may be con- 
sidered as a preference of the bump model as compared 
with the ramp model. 

Generally, the bump model can be considered as a 
more subtle and sophisticated version of the ramp model, 



13 



the former suggesting essentially the same predictions as 
the latter. Obviously, the ramp model is purely em- 
piric, while the bump model describes some dynamical 
processes in formation of vortices. Indeed, Fig. Off) 
presents two ramps, one is what is expected, and an- 
other one (dashed-dotted line) corresponds to a statisti- 
cally unlikely ramp. This statistical preference is so far 
not understood. It was suggested in Q] that the situ- 
ation is analogous to a shock wave formation: a shock 
would never develop in the shape like depicted with a 
dashed-dotted line, and it evolves only into "real ramp" , 
(solid line). This analogy is of little help, however, the 
typical turbulence being incompressible, and the shocks 
are therefore irrelevant. On the other hand, the bump 
should appear on descending part of the curve, where 
the derivative is negative, as depicted in Fig. 0(d), 
solid line, as opposed to "unrealistic" bump depicted 
with dashed-dotted line. This statistical preference can 
be explained, as opposed to the ramp- model statistics. 
Indeed, the ascending part of the curve (with positive 
derivative) corresponds to a structure function with vec- 
tor r parallel to the vertical direction, where the deriva- 
tive (= 2a) is indeed positive, while the descending part 
corresponds to r in horizontal plane (where the derivative 
is = —a < 0). Now, the bump corresponds to the vor- 
tex, and, of course, the corresponding velocity is suited in 
the horizontal plane. Therefore, the bump would appear 
only on descending part of the curve. 

Finally, let us discuss the quantitativ e est imates for 
the asymmetry. As mentioned in Sec. II A, the non- 



intermittency scenario does not account for the asym- 
metry. The latter is easily explained in extreme inter- 
mittency approach, see Sec. 

expression (^7|) does provide a real istic estimate of the 
skewness, 



The corresponding 
stimate of the 
However, this 



II B 



as noted at the end of 
estimate corresponds to extreme intermittency scenario, 
i.e., to unrealistic intermittency. Still, the asymmetry 
can be explained if we suggest that the vortex formation 
proceeds all the time. That is to say that the transi- 
tional stage mentioned in previous subsection is always 
present. If there is enough time for the perturbation to 
reach the Taylor microscale, then the vortex is complete, 
and therefore the intermittency is generated. But, con- 
sidering the fact that the intermittency is not that pro- 
nounced in fully developed turbulence, we conclude that 
the final formation of a vortex is a relatively rare event. 
More often, the vorticity of a perturbation is enhanced 
for a while, and then the generating cell ceased to ex- 
ist, and thus the generation stops. However, during this 
amplification, the vorticity of the generated vortex eas- 
ily exceeds the causing this generation strain, and then, 
according to (|26|), the third order structure function be- 
comes negative. This is valid for any particular scale 
r, if the vortex reaches it (and its vorticity exceeds the 
strain). This may account for the Kolmogorov law. 



VI. CONCLUSION 

It is known that production of vorticity, or vortex 
stretching, does not necessary implies that intermittency 
is generated. Indeed, in classical Kolmogorov picture 
K41, the energy cascade from large eddies to small actu- 
ally implies that small scale vortices are generated: be- 
cause the vorticity production proceeds homogeneously. 
In other words, each eddy decays into smaller eddies, each 
of which in turn decays into smaller eddies, and so on. 

As we argued above, the vortex stretching may be con- 
sidered as kinematic, in some limited sense, of course. 
For a passive scalar ip, the classical Kolmogorov picture 
is not much different from what is described in the pre- 
vious paragraph. That is, ip 2 is transported down the 
scales, so that ((V^) 2 ) is growing (which is similar to 
the vorticity production) . The energy cascade is gradual 
and continuous in this picture, meaning that the large 
scale fluctuations generate smaller scale perturbations, 
which in turn generate even smaller scales, and so on, 
until the perturbations reach diffusive scale where they 
eventually disappear. We can imagine a direct cascade of 
energy, without passing the intermediate scales in iner- 
tial range, directly from the large scales to the dissipation 
scale. To be more specific suppose, that the velocity tur- 
bulence is presented by large scale cells only (of the size 
£): like in Sec. |II Bj , but the velocity field is random in 
time, unlike time-independent, i.e., steady cells from Sec. 
[iTb|. Generally, maxima of S7ip will be generated in the 
vicinity of the stagnation lines and stagnation surfaces, 
in a time-scale ~ £/vg. The V?/> field will have the scale 
of Taylor microscale. However, if the life-time of the ed- 
dies T£ fa £/v£, as in regular turbulence, then, first, there 
would not be enough time to produce a sharp and con- 
centrated structure of the Vi/'-field during one turn-over 
time, and second, in the next life-time, the concentrated 
VV'-field structure is formed in a different place. As a 
result, the Vi/'-field is generated homogeneously in sta- 
tistical sense. In other words, there is no intermittency 
generated. Suppose now that the turbulent cells are per- 
sistent, or, at least, some of them are. Then, distinct 
VV'-field structures are generated, and they are as per- 
sistent as the cells themselves. The scenario described in 



Sec. [I B| of vortex generation is essentially the same. 

We presented an evidence of the presence of vortices 
in high Reynolds turbulence. This is clear from the fig- 
ures: the intermediate scale range (|32|), where these vor- 
tices are supposed to appear (see the beginning of Sec. 
Ill, and also |i|), is clearly noticeable. The main feature 
which distinguishes this range is a substantially enhanced 
intermittency. The vortices are aligned to the interme- 
diate eigen- vectors, ||, ||, [0. And, most important, 
these vortices account for negative skewness, that is, they 
result in a right asymmetry of the PDF. This picture pro- 
moted a development of the ramp model which is mod- 
ified now into the bump model. The latter not only ex- 
plains more features than the former, but also the bump 
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model seems to be dyn amic ally substantiated, unlike the 
ramp model, see Sec. 



VB 



The most important conclusion is that the intermit- 
tency is related to the asymmetry of statistical properties 
of turbulence, which was suggested in Q|. This give us 
a powerful tool for studying the intermittency: because 
the asymmetry is manifested in lower moments of struc- 
ture function, while the intermittency normally manifests 
itself in higher moments with poor statistics. 

Still, the whole picture is far from complete. As men- 
tioned in Sec. 



[II 



this simple scenario, outlined again 
above in this section, predicts too much intermittency. 
There are some suggestions in Sec. [II as to how to cure 



this difficulty. But, so far, the theory is unable to predict 
the numbers, e.g., we cannot predict the flatness from the 
asymmetry, or the other way. Nevertheless, we believe 
that the paper provides a next step in understanding the 
intermittency of turbulence. 
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